Phenomenological analysis of ATP dependence of motor protein 
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In this study, through phenomenological comparison of the velocity-force data of processive motor 
proteins, including conventional kinesin, cytoplasmic dynein and myosin V, we found that, the ratio 
between motor velocities of two different ATP concentrations is almost invariant for any substall, 
superstall or negative external loads. Therefore, the velocity of motor can be well approximated 
by a Michaelis-Menten like formula V = [ATP]fc(F)L/([ATP] + K M ), with L the step size, and 
k(F) the external load F dependent rate of one mechanochemical cycle of motor motion in satu- 
rated ATP solution. The difference of Michaelis-Menten constant Km for substall, superstall and 
negative external load indicates, the ATP molecule affinity of motor head for these three cases are 
different, though the expression of k(F) as a function of F might be unchanged for any external load 
F. Verifications of this Michaelis-Menten like formula has also been done by fitting to the recent 
experimental data. 



The processive motor proteins, including kinesin, 
dynein and myosin are essential for biophysical function- 
ing of eukaryotic cells 0, H- Due to the development 
of experimental instrument Q] , much accurate exper- 
imental data have been obtained Both conven- 
tional kinesin and cytoplasmic dynein move hand-over- 
hand along microtubules with step size 8.2 nm by con- 
verting chemical energy stored in ATP molecules into me- 
chanical works (ic3. [hi- 17|. Myosin (V or VI) also moves 



hand-over-hand but along actin filament with step size 
around 36 nm 0, [lsli^ . So far, there are many bio- 
physical models to understand the mechanism of motor 
proteins, including the flashing ratchet model (Tl], [21 , 22 1. 
Fokker-Planck equation 
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Meanwhile, more de- 



tailed mechanochemical models have also been designed 
to explain the experimental data, and get meaningful bio- 



ferent constant Km for substall, superstall and negative 
external loads. So the limit velocity in saturated ATP 
solution is V* — k(F)L, and the velocity corresponding 
to ATP concentration [ATP] can be obtained simply by 
multiplying V* by a constant [ATP] /([ATP] +K M ). 

For the sake of comparison, the velocity-force data of 
kinesin, dynein and myosin are plotted in Figs. [TJ [2] 
and[3fa). In Fig. QJa), the thick dashed line V\ is the 
velocity-force data of kinesin for [ATP] =1 mM obtained 
by Nishiyama et al [6], and the solid line V2 is for [ATP] 
=10 fj,M. One can easily see that there is only little differ- 
ence between the lines V2 and Vi/3.3. Similar phenomena 
can also be found for the velocity-force data of dynein 
and myosin obtained in 0, 10, [33], see Figs. [IJb,c,d). 



chemical parameters 
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In this study, by phenomenological comparison of the 
velocity- force data of different ATP concentrations, we 
found that the velocity of processive motor proteins 
can be described by a Michaelis-Menten like formula 
V = [ATP]/c(F)L/([ATP] + K M ), but might with dif- 
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Generally, these phenomena also exits for negative and 
superstall force cases, but might with different ratio con- 
stant, see Figs. 2] and [3ja) for data of kinesin obtained 
in Refs. Q, Q) Is [ ■ For the kinesin data in Q , the ratio 
constant is about 2.6 for F < 0, about 5.6 for < F < 7 
pN, and about 2.3 for F > 7 pN [see Fig. [S^a)]. For 
the data in [7j, the ratio constant is about 29 for F < 0, 
and about 16 for < F < 5 pN [see Fig. fflb)]. But for 
the kinesin data measured in Q, the constant 3.6 works 
well for both substall and negative external load [see Fig. 
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Data for kinesin from Nishiyama 2002 



Data for dynein from Toba 2006 



Data for kinesin from Carter 2005 



Data for kinesin from Block 2003 
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Data for myosin from Uemura 2004 




2 4 
F(pN) 

Data for Myosin from Tsygankov 2007 
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FIG. 1: For the positive substall external load cases, the ve- 
locity Vi of motor protein at low ATP concentration can be 
well approximated by the velocity Vi at high ATP concentra- 
tion divided by a constant, (a) For the experimental data of 
kinesin measured in j| i the velocity V2 of [ATP] =10 fiM can 
be approximated by Vi/3.3 with Vi the velocity of [ATP] =1 
mM. (b) For the data of dynein measured in [101 ]. velocity V2 
of [ATP] =10 fj,M can be well approximated by Vi/6.5 with 
Vi the velocity of [ATP] =1 mM. (c) For the data of myosin 
V measured in @|, velocity V 2 of [ATP] =10 can be well 
approximated by Vi/4.5 with Vi the velocity of [ATP] =1 
mM. (d) For the data of myosin V used in [32|] (derived from 
Q]), velocity V2 of [ATP] =1 /iM can be well approximated 
by 1/i/13 with Vi the velocity of [ATP] =2 mM. 



From the above observations about the experimen- 
tal data plotted in Figs. [1] and [2J one can see that 
the velocity-force relation of motor proteins satisfies 
V(F, [ATP]) = /([ATP])V;(F). Where V* = V*(F) 
is the velocity-force relation at saturated ATP concen- 
tration, and obviously V, can be written as V*{F) = 
k{F)L with L the step size of motor protein, and k(F) 
the force dependent rate to complete one ATP hydrol- 
ysis cycle (coupled with one mechanical cycle). The 
function /([ATP]) increases with [ATP], /(0) = and 
/([ATP]) = 1 with [ATP] -+ oo. A reasonable form of 
/([ATP]) is /([ATP]) = [ATP]/([ATP] + K M ) with a 
parameter Km which we called Michaelis-Menten con- 
stant. Finally, the velocity formula can be written as 
V(F, [ATP]) = [ATP]fc(F)X/([ATP] + K M ). 

To verify the above velocity-force formula, the force 
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FIG. 2: For general external load cases, the velocity V2 of 
kinesin at low ATP concentration can be well approximated 
by the velocity Vi at high ATP concentration divided by a 
constant, (a) For the data in [jj, the velocity Vi of [ATP] 
= 10 /iM can be well approximated by velocity Vi for [ATP] =1 
mM divided by a constant Const with Const=2.6 for F < 0, 
Const=5.6 for < F < 7 pN, and Const=2.3 for F > 7 
pN. (b) For the data in 0], the velocity V 2 of [ATP] =4.2 
fiM can be well approximated by velocity Vi for [ATP] =1.6 
mM divided by a constant Const with Const=29 for F < 0, 
Const=16 for F > 0. 

dependent expression of rate k(F) should be given firstly. 
Usually, the mechanical coupled cycle of ATP hydrolysis 
includes several internal states, here, as demonstrated 
in the previous mechanochemical model 27|, we assume 
that, in each cycle, there are two internal states, denoted 
by state 1 and state 2 respectively. 

one cycle 

..•5=*T^2FT *=*••• (i) 

Let Ui,Wi be the forward and backward transition rates 
at state i, then the steady state rate k(F) can be obtained 
as follows 27 



k= , (2) 

u\ + u 2 + wi + w 2 

The force dependence of rates Ui,Wi are assumed to be 

27] 



Ul(F) : 

W!(F) 



9+FL/k B T 



,jO e B-FL/k B T ^ 



u 2 (F)=u 2 e- e i FL ^ T , 
w 2 {F)^w° 2 e e - FL ' kBT . 



(3) 

Where ks is the Boltzmann constant, T is the absolute 
temperature, and Of' are load distribution factors which 
satisfy 8^+9^ + 9^+9^ = 1 . For this two-state model, 
one can easily get the following formula of motor velocity 
[ATP](uiu 2 — u>iw 2 )L 



V(F, [ATP]) 



' ([ATP] + K M ){ui +u 2 + wi+ w 2 ) ' 



(4) 
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Data for kinesin from Guydosh 2009 
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FIG. 3: (a) For the kinesin data measured in [J], the veloc- 
ity Vi of kinesin at low ATP concentration (10 fiM ) can be 
well approximated by the velocity Vi at high ATP concentra- 
tion (2 mM) divided by a constant 3.6, which is the same for 
both substall and negative external load. (b)Experimental 
data for conventional kinesin measured in [9J and the the- 
oretical prediction using the Michaelis-Menten like formula 
V = [ATP]fc( J F)L/([ATP] + Km)- The ATP concentrations 
are corresponding to [ATP] =1 mM (dashed line and squares) 
and 10 fiM (solid line and dots) respectively. The model pa- 
rameter K M is 15.8 (M for F < 0, 39.2 /xM for < F < 7 
pN, and 11.9 /tM for F > 7 pN, others are listed in Tab. |T] 




FIG. 4: Experimental data for cytoplasmic dynein ob- 
tained in [2| and myosin obtained in [5J (see also for 
the method to get the present values), and the theoreti- 
cal prediction using the Michaelis-Menten like formula V = 
[ATP]fc(F)L/([ATP] + Km)- (a) The experimental data are 
for [ATP] =1 mM (dashed line and squares) and 10 fiM (solid 
line and dots). The model parameter Km = 60.3 fiM for 
< F < 7 pN. (b) The experimental data are for [ATP] =2 
mM (dashed line and squares) and 1 [iM (solid line and dots) . 
The model parameter Km = 14.7 /iM for < F < 3 pN. 

The fitting results of the above velocity-force formula 
to kinesin data measured in [9] are plotted in Fig. Elb). 
In which, the Michaelis-Menten constant Km = 15.8 /iM 
for F < 0, K M = 39.2 /iM for < F < 7 pN, and 
Km = 11-9 /iM for F > 7 pN, other parameter values 
are listed in Tab. [I] Meanwhile, the fitting results to 
the dynein data measured in 10( and myosin data mea- 
sured in 0] are plotted in Fig. |4][a) and Fig. |4jb) (with 



TABLE I: Parameter values used in the theoretical predic- 
tions of the velocity-force relation for conventional kinesin, 
cytoplasmic dynein and myosin V: see Figs. EJb) and|lja)(b). 
The unit of rate u®,Wi is s _1 . 
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Michaelis-Menten constant K M = 60.3 fiM and 14.8 /iM) 
respectively, see also Tab. U for the corresponding pa- 
rameter values (34)]. Note, the step size used in the cal- 
culations is L — 8.2 nm for motor proteins kinesin and 
dynein, but L = 36 nm for myosin V. Certainly, the same 
fitting process can also be done to other experimental 
data. The plots in Figs, ^[b) and H] indicate that, the 
experimental data of motor proteins can be well repro- 
duced by the Michaelis-Menten like formula (j4]), so our 
phenomenological analysis about the ATP dependence of 
motor motion is reasonable. 

In summary, in this study, the ATP dependence of 
motor proteins is phenomenological discussed. Based 
on the recent experimental data and basic numerical 
calculations, we found the motor velocity can be well 
described by a Michaelis-Menten like formula V — 
[ATP]k{F)L/([ATP] + K M ) with force dependent rate 
k(F) at saturated ATP. The different values of Km for 
substall, superstall and negative external load imply, 
the ATP molecule affinity of motor head might be dif- 
ferent for these three cases, but the basic mechanism 
in each mechanochemical cycle (either forward or back- 
ward) might be the same. An obvious conclusion from 
our Michaelis-Menten like formula is that the stall force, 
under which the mean motor velocity is vanished, is in- 
dependent of ATP concentration 0, (3 ■ 
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